function [yx, yz, xz, profile] = plot_V(V_1)
    deta_z = 125/20/15;
    deta_x = deta_z;
    % V_1 = imresize3(V_1, [size(V_1, 1), size(V_1, 2) floor(size(V_1, 3) /( 125/40/15))]);
    % V_1 = imresize3(V_1, [size(V_1, 1), size(V_1, 2) length([-100:deta_z:100])], 'method', 'linear');
    figure,
    yx = squeeze(max(V_1, [], 3));
    
    yz = squeeze(max(V_1, [], 2));
    xz = squeeze(max(V_1, [], 1));
    imagesc(yx), axis off, axis equal, title('yx');
    
    figure,
    subplot(1, 2, 1), imagesc(yz), axis off, axis equal, title('yz');
    subplot(1, 2, 2), imagesc(xz), axis off, axis equal, title('xz');
    % figure,
    % imagesc(squeeze(max(V_1, [], 2))), axis off, axis equal, title('zy');
    %%
    % figure,
    % % subplot(1, 2, 1), plot(squeeze(V_1(28*15+8, :, floor((31+15)/61 * length([-30:deta_z:30]))))), title('profile of x');
    % subplot(1,2,1), plot([-100:deta_z:100], squeeze(V_1(36*15+8, 38*15+8, :))), title('z profile of single beads z= 10');
    % subplot(1, 2, 2), plot([-100:deta_z:100], squeeze(V_1(28*15+8, 38*15+8, :))), title('z profile of 3 beads z = [6, 10, 15]');
    % 
    % profile(1, :) = [-100:deta_z:100];
    % profile(2,:) = squeeze(V_1(28*15+8, 32*15+8, :));
    % y =  squeeze(V_1(36*15+8, 32*15+8, :));
    % x = [1:size(V_1, 3)]*deta_x;
    % [fitresult, gof] = createFit(x, y);
    % fitresult
    
    % 2.355*std = FWHM
end